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SUMMARY 


A number of accepted techniques exist to invert radiant transfer data. 

Each technique involves certain assumptions and has inherent limitations. 
Derived here is a new iterative technique which recovers density profiles in a 
nonhomogeneous absorbing atmosphere. The technique is based on the concept of 
factoring a function of the density profile into the product of a known term 
and a term which is not known, but whose power series expansion can be found. 
This series converges rapidly under a wide range of conditions. 

A demonstration example of simulated data from a high-resolution infrared 
heterodyne instrument is inverted. For the examples studied, the technique is 
shown to be capable of extracting features of ozone profiles in the troposphere 
and to be particularly stable. 


INTRODUCTION 

One of the fundamental problems in the theory of atmospheric remote sensing 
can be posed as follows; If the radiation emitted from an object whose emissive 
characteristics are known is measured in each of a number of instrument channels 
after the radiation has passed through an absorbing atmosphere with known or 
assumed optical properties, what can be inferred from these measurements about 
the vertical distribution of the absorbing species in the atmosphere? For exam- 
ple, one may want to infer the concentration profiles of an atmospheric gas or 
the temperature profile of the atmosphere from measurements made by an instru- 
ment looking toward the Earth from a satellite and sensing radiation emitted 
from the Earth in several frequency bands or at several different viewing 
angles . 

Considerable work has been done in the field of inversion theory (refs. 1 
to 5) . As a result, a number of reliable techniques are available to "invert" 
radiant transfer data. The purpose of this paper is to derive a new iterative 
technique for inferring concentration profiles from such a set of measurements. 
In order to examine some of the convergence properties of this method, it is 
used to invert two numerical demonstration examples. The results are compared 
with those of other standard inversion techniques. 

The present method is very similar to that developed by W. L. Smith 
(ref. 1), the difference being in the mathematical form of the "weighting 
functions" described herein. Although neither method is restricted to a 
homogeneous atmosphere, both methods reduce to the same algebraic form for 
the special condition of nadir viewing in a homogeneous medium. 



SYMBOLS 


B Planck function 

f general function in equation (1) 

g,h general functional form 

I radiation spectral intensity 

J total number of instrument channels 

K absorption cross section 

n n( 2 ), number density 

p pressure at altitude z 

Pq sea-level pressure 

S weighting integral defined by equation (AT 8) 

T transmittance 

U optical mass 

W weighting function 

y In (p/Po) 

z altitude measured frcsn sea level 

z' variable of integration 

e function defined by equation (A6) 

y cosine of angle emerging ray makes with local vertical 

T optical depth 

()) instrument response function 

Subscripts; 

j channel number, 1, 2, . . J 

TOP top of atmosphere, or satellite altitude 

X wavelength 
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Superscripts: 

m iteration number 

meas measured value 

* particular altitude in interval [0,z^Qp] 


PRELIMINARY REMARKS 

In the present approach, the intensity at the top of the atmosphere 
IX(z^Op) is related to the concentration throughout the atmosphere by an 
equation of the form 

f^TOP r , 

i\(zTOP'^) = Bx (z^op) + ) dz* fxLKx (z') (z ' ) ,y,n (z ') J (1) 

in which the quantity to be determined is the vertical distribution of the 
absorbing species n(z). If equation (1) is put into an iterative form, the 
difference between I\ calculated using the (m + 1 ) th iteration for n(z) 
and that using the mth iteration is 


AIx(ztoP'M) = Ix^"''*'^^ (z^PfU) “ 

dz' {fx[Kx(z'),Bx(z'),M,n ("'+■') (z’)] 

~ fx[Kx (z ' ) ,Bx (z ' ) ,M,n"* (z ' ) ]^ (2) 

A procedure has been found for factoring the braced term so that it can be 
written as the product of two functions, 

g(z,vi) • h[z,e (z) ,y] 

where g(z,y) is a known function of n"*(z) and h[z,e(z),y] is an unknown 
function which depends on n ) (z) , through the term e(z). It is shown in 
the appendix that e(z) is bounded by a small parameter for all values of z. 
Thus, by expanding h[z,e(z),y] in a power series in £(z) and retaining only 
the linear terms, equation (2) can be put into the form 

f^TOP 

Aix(ztop»M) = dz' ex(z') Wx(z’) (3) 

where Wx(z) is known and is usually referred to as the "weighting function." 
In this method, the weighting function Wx(z) is given by equation (A21 ) , 


pzrpop 
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<3Bx(z) 


( 4 ) 



where is the optical depth, T\ is the transmittance, and B;^ is the 

Planck function. 

The weighting function concept is one which is used throughout the field 
of inversion theory. This function varies with altitude and essentially gives 
the user a picture of where the information is being obtained. If, for example, 
the weighting function peaks in the stratosphere at a certain frequency and 
angle, then practically all of the measured information is coming fron high 
altitude. If information in the troposphere is required, look for weighting 
functions which peak in the troposphere. As discussed in the references 
(refs. 1 to 5) , it is important to have a reasonably good distribution of 
weighting function peaks with altitude. To avoid ambiguous results, it is 
important that the peaks overlap as little as possible. 

As can be seen from equation (4), the weighting functions in the present 
method vary with y and X. The functional dependence on y not only occurs 
in the factor 1 /y but also in the transmission term 

Tx(z,y) = exp[-Tx(z)/y] (5) 

The functions xx(z), Tx(z,y) and Bx(z) all depend on the wavelength of the 
radiation being measured, since these terms all depend heavily on the molecular 
absorbing properties of the medium of transmission. 

Since the intensity at the top of the atmosphere is measured by an instru- 
ment, the quantities such as T and x being used in the weighting functions 
are more appropriately thought of as averaged quantities over the band-pass 
ranges of the instrument channels. Thus, for the jth channel, assuming a rec- 
tangular response function 


( 6 ) 


and so on for all wavelength dependent quantities. The width of the channel is 
represented by AX, with ^lower ^upper representing the lower and upper 

wavelengths of the channel and representing the instrument response func- 

tion. If a high resolution instrument is being used, Tj and Tx will be 
reasonably interchangeable for any wavelength in the channel. The narrower the 
channel, the more accurate the approximation will be. 


upper 


4)(X) Tx(z,y)dX 


Xj (z) = ^ J 
J AX 


4)(X) xx(z)dX 


lower 
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A RETRIEVAL PROCEDURE 


Measurements of the radiation emanating fran and passing through the atmo- 
sjAiere are made in several discrete channels (j = 1 , 2, . . . , J) , each with its 
own finite field of view, bandwidth QAXj = (^j) upper “ (^j)lowerD' viewing 

angle with respect to the outward normal (cos“^ y) . The problem faced by the 
researcher is determining how to infer, from these measurements, a vertical 
profile of a trace gas which is optically active to all channels and which has 
a vertical distribution consistent with all of the measurements to within the 
measurement error. The inversion scheme presented here provides a means for 
accomplishing this. 

The inversion scheme outlined here is derived in detail in the appendix. 

The procedure, suitable for use on a high-speed computer, can be reduced to the 
following sequence of steps: 

1. Assume an initial gas profile n®(z) as a first estimate 

2. In each channel, calculate the upwelling radiance at the top of the 

atmosfiiere by using equation (A2) , 





'^Ztop ' ) 

i dz' 

'O ‘^z* 



^^TOP 

X exd 

1 1 

1 K-j(z") n"»(z")dz" 


_■ yj - 

'o 


(7) 


3. Calculate the weighting function for each channel by use of 
equation (A21 ) , 

1 dBj (z) 

W-i"'(z) = - T^m(z) T-i"*(z) 

^ y J J dz 


( 8 ) 


where ij™ is the optical depth predicted from the measurements in 
the jth channel on the mth iteration and Tj™ is the corresponding 
transmittance computed from this value of the optical depth. 


4. For each data channel, calculate a new gas profile n j ) (z) 

of equation (A20) , 


by use 


nj (!”+■') (z) 


„ / ''ij" 

n™(z) 1 1 + 

\ . 


(9) 
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where 


Alj"> = [lj"*®^®(ZTOP.Uj) - Ij">(ZTOPrVij)] 


no) 


and Sj™ is computed from equation (A18) 


1 r^TOP dBj(z') 

= — \ dz> T^"‘(z•) T4"1 (z') 

^ Uj -'o dz' ^ ^ 


(11) 


5- The final estimate of n(z) for each altitude is obtained from a 

weighted average of the individual profiles, given by equation (A22) , 


max 


n ) (z) = 


Vj. (m+1) ( 2 ) n-j ) (z) 

j=1 

Jmax 

> Wj ) (z) 

j=l 


( 12 ) 


There are a number of points worth discussing in connection with this 
algorithm. If the weighting functions are well separated with little or no 
overlap (which seldom happens in practice), then the best estimate for n(z) 
in a given altitude region is that given by the channel with a nonzero weight- 
ing function in that region; that is, if the jth channel is the only one which 
has a nonzero weighting function at sane altitude z* , then equation (A22) 
reduces essentially to 

0 + 0 + . , . + nj (z*) • Wj"*(z) + 0 + . . . 

n(z*) = — ^ = nj(z*) (13) 

0 + 0 + . . . + Wj"*(z*) + 0 + . . . 

Thus, in the altitude regions where there is some overlap of the weighting 
functions, the weighting procedure provides a reasonable scheme to form a physi- 
cally meaningful average of the information from the contributing channels. 


In step 4 of the procedure, the step size 


Alj"*(ZjQp,y j) 


is not bounded 


c . in 

^3 


mathanatically . It has been found through experience that the rate of conver- 
gence is frequently enhanced by restricting the step size to 0.5. This value 
is arbitrary, but it has been found through experience to be useful. Although 
limiting the step may require a few more iterations in some calculations, it 
was found expedient to approach the solution by more moderate size steps. In 
some calculations, limiting the step size may even prevent divergence. 
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It is worth repeating that this technique is closely related to that of 
W. L. Smith (ref. 1), as described by Seals in reference 6. The only major 
differences in the algorithms are as follows: 

Weighting function Wj"*(z) - Given in the present technique as 


1 dBj (z) 

W-i(z) = - T^n»(z) T-i^Cz) (14) 

J y ^ J dz 


it is given in the Smith technique (using present notation) by 

dBj (z) 

W^">(Z) = u"'(z) K-;(z) T-;"'(Z) (15) 

^ ^ dz 

where U"*(z) is the optical mass path given by the mth iteration, or by Seals 
as 

9Bj (y) 3Tj"'(y) 

Wj">(y) = — (16) 

9y 9y 


with y = In (p/Pq) when p is the pressure at level y and Pq is the 
pressure at sea level. 

S integral - Given in the present technique by 
1 P^TOP dBj(z) 

Si = - \ dz Tim(z) Ti"»(z) (17) 

^ y Jq ^ ^ dz 

it was given in reference 6 by 


S 


j 


i 


U">(y) 


ey 


9Bj (y) 

9y 


9Tj">(y) 

9y 


dy 


( 18 ) 


Thus, if the Smith algorithm is already available in program form, the present 
algorithm can very easily be added to the program, giving the user an additional 
option for inverting data. 


NUMERICAL EXAMPLE 

The inversion simulation was performed in three major steps: 

1 . A study was made to find a set of narrow data channels which was appro- 
priate for a high-resolution infrared heterodyne spectrometer and was reasonably 
free of the effects of interfering gases. 
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2. A set of synthetic "measurement data" was created using an upwelling 
radiance program which evaluated equation (Al ) for an assumed nominal density 
profile of the trace gas, ozone. 

3. The algorithm described in the section "A Retrieval Procedure" was then 
used to invert the measurement data in an attempt to recover the density pro- 
file. No measurement error was assumed in these studies. 

Detailed descriptions of these steps follow. 


Channel Selection 

To find appropriate wavelengths and angles at which to locate the instru- 
ment channels, the 1966 U.S. Standard Atmosphere (ref. 7) and a simple absorp- 
tion model were used in a study to calculate the total transmission through a 
homogeneous slab where pressure, temperature, and thickness are specified. 

This computer program, developed by R. K. Seals, Jr., of the NASA Langley 
Research Center and L. L. Gordley of Systems & Applied Sciences Corporation, is 
unpublished. It essentially calculates the monochromatic transmissivity at a 
number of discrete frequencies and is based on the assumption of one basic 
absorbing species and up to five interfering gases. Line parameters from an 
updated version of the McClatchey tape (ref. 8) are used. The output consists 
of two curves of transmittance versus wavelength - one for the basic absorber 
alone and one for the basic absorber plus the interfering gases. An example of 
this output is shown in figure 1 for the wave-number range of 1041.5 cm~^ to 
1044.0 cm~^ and in figure 2 for the range of 1050 cm“^ to 1060 cm”^ . In both 
figures, the basic absorbing gas was ozone and the interfering gases were 
carbon dioxide, water vapor, and methane. (See table I.) The water vapor 
continuum was not included in this model. 

Several regions of minimum interference can easily be seen in these fig- 
ures. To take advantage of these regions, the five channels shown in table II 
were chosen for the simulated soundings. 


Synthetic Data 

An upwelling-radiance program was developed to solve equation (Al ) for a 
large number of monochromatic frequencies. The results were then integrated 
over the finite bandwidth of each channel to produce the "measurements." A 
rectangular instrument response function was assumed. The neutral composition 
of the atmosphere used is presented in table III. 

In evaluating equation (Al ) , a plane-parallel atmosphere was used. This 
was broken into 50 layers, each of 2-km thickness. The Lorentz profile was 
used to compute the absorption cross section below an altitude of 34 km, and 
the Voigt profile, calculated from the computational algorithm of reference 9, 
was assumed above this altitude. 
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Inversion 


With the atmospheric upwelling-radiance program, all pertinent physical 
variables, such as optical depth and transmissivity, can be calculated sepa- 
rately for each species at each altitude and wavelength. After ozone was 
chosen for the gas to be profiled, several instrument channels were selected. 

The weighting functions for these channels were computed by equation (A21) and 
are pictured in figure 3. Note that peaks occur below the ozone bulge at 21 km 
(fig. 3) . The dependence on 1/y is evident in channels 3 and 5 where the 
same wavelength was chosen, but different angles were used. Channel 3 peaks at 
a slightly higher altitude, since its ratio of optical depth to M is greater 
than that of channel 5 at all nonzero optical depths. 

A number of profile retrievals were made corresponding to different initial 
guesses of the ozone profile. Two estimated profiles are shown in figure 4. 
After some experience was gained, it became apparent that an interesting feature 
of the technique presented herein was its stability, particularly in comparison 
with the closely related formulation of the technique of Smith in reference 6. 
However, it must be emphasized that, for the application of the present tech- 
nique, the channel selection was made to spread the peaks of the weighting 
functions (eq. (A21 ) ) as evenly as possible over the altitude range. These 
channels would not necessarily be suitable for the application of the Smith 
technique which used the weighting functions of reference 6. These weighting 
functions, for the same set of conditions as stated above, are shown in fig- 
ure 5. These all peak above the troposjAiere for the stated conditions and 
hence all provided "information" from the same region of the atmosphere. It is 
possible that a more nearly optimum selection of channels tailored to the Smith 
weighting function would yield results similar to those obtained by the present 
method. The following comparison of numerical results should, therefore, be 
viewed as an application of both methods to the same instrument and not as a 
comparison of optimum retrieval techniques. 

The three ozone profiles used to simulate the radiance measurenents for 
each of the five channels are shown in figure 4. The curve labeled "First 
Guess for Case 1" is the initial profile used in both the present and the Smith 
retrieval methods. The results of the retrieved profiles are shown in figure 6. 
The solid line shows the true profiles, the dashed lines show the results of 
the present method after 30 and 54 iterations, while the dotted line shows the 
results of the Smith technique after 63 iterations. The present method recovers 
the profile below the tropopause reasonably well and does not fare too badly in 
recovering the bulge and the trailing profile above the tropopause. There was 
no substantial improvement after 30 iterations, possibly due to accumulation of 
round-off error (ref. 2) . Using the same data, the Smith technique failed to 
achieve low residual values in its iterations. The residual error in this sim- 
ulation is defined to be 


Residual™ 



channels 


(Measured Radiancej 


Calculated Radiancej™) 


2 


The residuals from both methods are shown in figure 7. Those from the 
present method varied quite smoothly and stabilized at about 41 iterations. 
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The Smith residuals also seemed to stabilize somewhat at this point, at a level 
about an order of magnitude higher than those of the present method. The 
initial rate of convergence is about the same for both methods up to about 
15 iterations. 

It is worth emp^iasizing that, for initial guesses this far from the true 
value, the limitation of the maximum step size to |Ai/s| =0.5 was mandatory. 
Neither the Smith technique nor the present method converged without the use 
of this criterion. Negative number densities and extremely large radiances 
resulted when the step size was allcwed to vary without bound. 

The results just discussed were for a retrieval where the initial profile 
underestimated the true atmospheric profile at all altitudes. A profile 
retrieval using the present method is shown in figure 8. The initial profile 
has the same shape as the true profile, but everywhere overestimates the pro- 
file as shown in figure 4. The retrieved results are shown in figure 8 for the 
38th iteration, and the residuals are shown in figure 9. The step size inhibi- 
tion was used. Once again the present results reproduced the true profile 
(fig. 8) reasonably well, particularly below the tropopause. The initial 
behavior of the residuals was more erratic, but they eventually smoothed out 
and reached a value slightly higher than that of the initial profile. 


CONCLUDING REMARKS 

In this study a new iterative technique is derived which is capable of 
efficiently inferring density profiles fran radiant intensity measurements made 
in a nonhomogeneous absorbing atmosphere at an arbitrary angle with the outward 
normal. Stability has been demonstrated for the examples presented, and expe- 
rience indicates that the technique converges over a wide range of conditions. 
Instrument channels were selected to take advantage of the characteristics of 
high-resolution infrared heterodyne spectrcmeters which are capable of extract- 
ing information about ozone profiles well below the 21-km density profile bulge. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 19, 1979 
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APPENDIX 


A DERIVATION OP THE RETRIEVAL ALGORITHM 

In the present derivation, the assumption is made that the tenperature 
distribution through the atmosphere, as well as the density distribution of 
all optically active species other than the one under investigation, is known. 
Under conditions of thermodynamic equilibrium for a nonscattering planar atmo- 
sE^iere, the radiant transfer equation for a single absorbing species relates 
the intensity of radiation in channel j to altitude z as follows; 


Ij (ztop»V* j) 


Bj(0,Hj) 


exrt - 



Kj (z') 


n (z ' )dz 


1 P^itop 

+ — \ dz ' Bj (z ' ) 

Mj ^0 


X Kj (z ' ) n (z ' ) 



1 P^TOP 

— K.s(z") n(z")dz" 

Uj 


(Al) 


where the only sources of radiation are the Earth (which is assimed to radiate 
as a black body) and the atmosphere. 

Integration of the second term of equation (Al ) by parts yields 


I j (Z top»M j) 


Bj (z^oP'V* j ) 


pz,pOP dBj(z') 


1 p^^TOP 

1 dz' 

ejq) 

\ K-i(z'') n(z")dz" 

Jq dz' 


L "j J 


(A2) 

If equation (A2) is considered to be the equation governing the intensity 
of the mth step of an iterative solution for n(z), the difference between the 
intensity functions for the mth and (m + l)th iteration is 


pz^op dBj(z') 

j.(m+l) _ j.m = - \ dz' iexpi 

^ ^ dz' 


1 p^TOP , , 

— \ K-;(z") n (z")dz" 


- exp 


1 p^^TOP 

pj 4' 


Kj(z") n">(z")dz" 


(A3) 
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APPENDIX 


Define the optical depth Tj (z) at altitude z by the relation 
pZrpop 

Tjin(z) =J Kj(z') n»“( 2 ')dz' (A4) 

and substitute this into equation (A3) . As usual in iterative methods, the 
assumption is made that number (the radiation at the top of the atmo- 

sjiiere in the jth channel of the iteration following the mth) can be replaced 
by the value of the intensity in the jth channel as measured by the instrument. 
Thus, 


pz^op 

Aljin H IjineaS(z^QP^jj ) _ i.m(z^Qp^p ) = ^2' 

t/n 


dBj(z') 
dz ' 


X 



(A5) 


Define the function €j"'(z} such that 
Tj (m+l) (z) 

= 1 + £j"‘(z) (A6) 

Tj"*(z) 

Then the expression in parentheses in equation (A5) can be simplified to 


exp|^-Tjn>(z')[l + ejro(z')^^ljjJ - exp[-T jn>(z ' )/)j j] 


= exp^-Tjin(z')]/M^ ^exp^[-£j">(z') Tj™(z')] 


Uj \ - 1 


(A7) 


After completion of the mth iteration, equation (A7) reduces to the product 


of two terms. The first term, exp \ [|-Tj’"(z ’) 


^h} ■ 


is completely determined. 


The second term, containing the parameter (z), is unknown, but its power 

series expansion will converge rapidly if [^Tj™ (z ' ) j is bounded and if 
Gj™(z') can be shown to be small enough throughout 0 ^ z' < ZipQp such that 
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APPENDIX 




« 


(A8) 


It has been found through experience (e.g., in the numerical demonstrations 
presented in the body of this paper) that Gj"*(z) is of the order of 1 0~^ or 
10“^ and becomes even smaller as Al™ ->■ 0. Thus, since the product of et 
must be small, optical depths in excess of 10 are routinely within the range of 
validity of this approach. 


By making the power series expansion and retaining linear terms, the right- 
hand side of equation (A7) can be written as; 


Use of equation (A9) results in equation (A5) reducing to 


(A9) 


1 f'^TOP dBj(z') 

Ai-m = — \ (321 T-;n'(z') T-!"*(z') 

^ yj '-'o dz' 


(A10) 


where 


Tj">(z) = exp[-Tj"'(z)/yj] 


In order to show that Gj^'Cz) is indeed "small," and to justify further 
manipulation of equation (A10), the central assumption is made that the Xj™(z) 
converge uniformly on the interval 0 ^ z < Zijop as m ->■ oo. With this assump- 
tion, the algorithm can be derived. 

That Ej^**(z) is arbitrarily small on the interval 0 ^ z < ZrpQp follows 
frcm rewriting equation (A6) as 


•Y . (m+1 ) (z) - T^™(z) 

^ = Gjn'(z) 

Tj"‘(z) 

Since Tj™(z) converges uniformly on 
exists an integer N such that for all m 
contained in the interval Qo,ZrpQp], 


(All) 

Qo,Zijiqpj, this implies that there 
larger than N and for all z 


]Tj (z) - Tj"»(z) I < 6* 
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APPENDIX 

for a 6' still to be determined. Since z < z.pQp, Tj™(z) > 0. Consider any 
Y > 0. By defining 6' = Y'Tj"*(z*) where Xj™(z*) is the maximum value of 

Tj™(z) on [|0,zipop3' then £j™(z) < y for any y provided N is sufficiently 
large . 


The uniform convergence of Tj”*(z) will also be used to show that equa- 
tion (A1 0) can be rewritten (with an arbitrarily small error as m ^ oo) as 
follows: 


Al 


1 dBj(z') 


•m.L r 

3 u. J. 


dz ' 




dz' 


e^"'(z') T^">(z’) T^™(z’) 


ej"'(z) pz^op dBj(z') 


Mj *^0 


dz ' 


dz 


Tji"(z') Tj"'(z') + A(z) 


(A12) 


where A(z) is the error created by bringing e(z') outside the integral and 
will be shown to be arbitrarily small for all values of z in f0,z.j.Qp^ (and 
can be neglected for computational purposes) . 


Consider the difference 


A(z) = — 


1 p^TOP j ) 


J dz' j Tj*"(z') Tj™(z') [|£j"*(z') - £j"’(z)J 


(A13) 


where )Jj ^ 0, dBj(z')/dz', 'Tj"'(z'), Tj™(z') are all physical quantities 

he 


and are hence bounded so that 


1 

P-i 


dB^ 


dz 


TT 


< M 


0 ^ z < z 


TOP 


for some M. Thus, 


i^zTOP 

A(z) < dz'M[_£jn>(z ') - £j">(z)3 


(Al 4) 


But since it has just been shown that for any z in the interval 


[^0 , z^opj I 


£j"*(z) < Y for arbitrary y > 0 

if m is sufficiently large. For ejn'(z'), choose Y = S/Mz^op- For Gjn>(z), 
choose Y = 5/2MzrpQp. Then, 
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C^TOP 6 

A (z) < \ dz'M 


■'0 2MZq;QP 2 

Since 5 is arbitrary, equation (AlO) can be rewritten as 


(AT 5) 


ej™{z) pzTOP dBj(z') 

Ai-m = I ^2' — T-;™(z') 

... dz' ^ ^ 


Vij *^0 

with negligible error 
For any z, then 


where 


Alj"* = 


Tj (z) - Tj">(z) 


Tj™(z) 


Q t m 


(A16) 


(AT 7) 


1 f^TOP ‘^Bj (z ' ) 

S.m = — I j 32 ' T^">(z') T^">(z') 

J II.1 Jn -J 




(A1 8) 


Thus, 


I 

T- j (m +1 ) (2) = Tj"'(z)jl + 


c .m 
^1 


(AT 9) 


By differentiating equation (AT 9) with respect to z and using equation (A4) 
to eliminate the derivative, this can be put into the form of an iterative 
algorithm for the number density n(z). 


I 

nj (m+1 ) (z) = nj^Cz)!! + 


(A20) 


c « m 
^3 


The integrand of Sj in equation (Al 8) is the weighting function 


Wj">(z) = — Tj">(z) Tj™{z,y) 


dBj (z) 
dz 


(A21) 
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APPENDIX 


With these weighting functions, the average n(z) can be defined as 


max 


n ) (z) = 


> Wj ) (z) n j ) (z) 

j=l 

3 max 

> Wj ) (z) 

j=l 


(A22) 


In the iterative algorithm described in the text, the average n(z) will be 
used in equation (A20) rather than nj (z) . This represents a physically 
reasonable interpolation scheme between altitudes with weighting function 
peaks. It should be noted that since Alj>" and Sj™ are height independent, 
equation (A20) acts to alter the magnitude of n(z), but not its profile shape. 
Equation (A22) does all the profile shape alteration. 
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TABLE I.- HOMOGENEOUS SLAB USED AS BASIS FOR SELECTING ABSORPTICW CHANNELS 


Path length 70 km 

Pressure 50 kPa 


Tanperature 


240 K 


Composition; 


O 3 0.1140 ng/cm^ 

CO 2 97.20 ng/cm^ 

H 2 O 120.0 ng/cm^ 

NH 3 0.5720 pg/am^ 


TABLE II .- CHANNELS SELECTED FOR SAMPLE CALCULATIONS 


Channel, j 

Initial wavelength, cm ^ 

Final wavelength, cm ^ 


1 

1 042.20 

1 042.23 

1 

2 

0 

00 

1 041 .90 

1 

3 

1041.58 

1041 .61 

.707 

4 

1 050.85 

1050.88 

.500 

5 

1 041 .58 

1041 .61 

.500 


1 8 



TABLE III 


(a) Ph 


Altitude, 


km 


Temperature, 


K 


99 

1 85.86 

95 

176.59 

91 

167.32 

87 

165.06 

83 

168 .70 

79 

1 78.50 

75 

1 96 . 1 0 

71 

21 3.70 

67 

231 .70 

63 

246.90 

59 

259.54 

55 

269.30 

51 

274.42 

47 

273.1 5 

43 

264.95 

39 

255.05 

35 

245.15 

31 

235.75 

27 

227.50 

23 

222.80 

1 9 

21 8.00 

1 5 

21 5.70 

11 

228.90 

7 

254.70 

3 

279.20 

1 

289.70 


OF ATMDSPHERE 


Pressure, Pa Density, kg/m^ 


0.0361 
.0706 
.138 
.299 
.667 
1.45 
2.99 
5.70 
1 0.39 
18.20 
30.94 
51 .53 
84.17 
1 37.26 
226.24 
380.18 
652.55 

1 145.57 

2 054.87 

3 753.08 
6 947.86 

13 030.40 
24 267.30 
42 047.70 
71 069.40 
90 240.00 


0.601 

X 

10-6 

1 .26 

X 

10-6 

2.63 

X 

10-6 

21 .9 

X 

10-6 

45.7 

X 

10-6 

23.2 

X 

10-6 

43.3 

X 

10-6 

76.0 

X 

10-6 

.128 

X 

10-3 

.21 2 

X 

1 0-3 

.345 

X 

10-3 

.560 

X 

10-3 

.909 

X 

10-3 

1.50 

X 

1 0-3 

2.60 

X 

1 0-3 

4.63 

X 

10-3 

8.77 

X 

10-3 

1 5.7 

X 

10-3 

29.2 

X 

10-3 

55.0 

X 

10-3 


.104 

.195 

.365 

.590 

.909 

1.112 




table ui.~ 


Concluded 


<b) »«nber density of 


absorbing 


™olecules/m3 

4-519 X TO? 3 
8-228 X 1 q1 3 
^ - 497 X 1 q1 4 
2-883 X 1 qT 4 
5-654 X jq14 
^ - TOO X 1 q1 5 
2-096 X 1 q1 5 
3*992 X to15 
8-248 X 7o 15 
^-757 X 70^6 
3-773 X igl 6 

8.061 X 1 q1 6 

^ • 740 X 7 0^ 7 
3-754 X tq17 
8.093 X 7o77 
7 - 787 X 7 q7 8 
^* 7 74 X 7 q7 8 
9-780 X 7q78 

^ -293 X 7q7 9 

T .412 X 7q7 9 
2- 478 X 7 q7 9 
7 -797 X 7q20 I 

7-304 X 7o27 1 

7 - 078 X 7 q22 j 

9-309 X 7q22 / 

2-735 X 7 o23 I 


speci es 


3/ molecules/m3 

8-623 X 7 q7 ] 
2.104 X 1 q1 2 
5- 7 36 X 7 q7 2 
^ • 336 X 1 q1 3 
3.740 X 7q73 
^ • 035 X 7 q 7 4 
2*959 X 7 o7 4 
3*459 X 7 q74 
2-029 X 1 q1 5 

^•597 X 7q75 
^• 027 X 7 q76 
2-762 X 1 q1 6 
4.576 X 7 q7 6 
^*085 X 7q77 
2.674 X 7 q7 7 
8* 764 X 7 q7 7 
7 - 7 62 X 7q7 8 

2- 7 92 X 7 q7 8 

3- 073 X 7q78 
3- 726 X 7 q7 8 
3-827 X 7 q7 8 
2- 257 X 7 q78 

^ * 324 X 7 q7 8 
^•025 X 1 q1 8 
3*707 X 7q77 j 

3.022 X 7q7 7 / 


1041.5 


1042,0 


1043.5 


1044 


1042 


5 

Wave number, on" 


1043.0 

■1 


Figure 1.- Transmittance as function of wave number for homogeneous slab from 1041.5 to 1044.0 cm 
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Figure 2.- Transmittance as function of wave number for homogeneous slab from 1050 to 1060 cm“^ . 
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Figure 



3.- Set of weighting functions for ozone determined by present method 





Altitude, 



Figure 5.- Set of weighting functions determined by Smith method 
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Figure 6.- Results of inversion algorithms at various iteration steps, first 

demonstration case. 
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Figure 7.- Residuals for first demonstration 


case 
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